Deciphering Bacterial and Archaeal Transcriptional Dark Matter and Its Architectural Complexity.

Transcripts are potential therapeutic targets, yet bacterial transcripts remain biological dark matter with uncharacterized biodiversity. We developed and applied an algorithm to predict transcripts for Escherichia coli K12 and E2348/69 strains (Bacteria:gamma-Proteobacteria) with newly generated ONT direct RNA sequencing data while predicting transcripts for Listeria monocytogenes strains Scott A and RO15 (Bacteria:Firmicute), Pseudomonas aeruginosa strains SG17M and NN2 strains (Bacteria:gamma-Proteobacteria), and Haloferax volcanii (Archaea:Halobacteria) using publicly available data. From >5 million E. coli K12 ONT direct RNA sequencing reads, 2,484 mRNAs are predicted and contain more than half of the predicted E. coli proteins. While the number of predicted transcripts varied by strain based on the amount of sequence data used for the predictions, across all strains examined, the average size of the predicted mRNAs is 1.6-1.7 kbp while the median size of the predicted bacterial 5'- and 3'- UTRs are 30-90 bp. Given the lack of bacterial and archaeal transcript annotation, most predictions are of novel transcripts, but we also predicted many previously characterized mRNAs and ncRNAs, including post-transcriptionally generated transcripts and small RNAs associated with pathogenesis in the E. coli E2348/69 LEE pathogenicity islands. We predicted small transcripts in the 100-200 bp range as well as >10 kbp transcripts for all strains, with the longest transcript for two of the seven strains being the nuo operon transcript, and for another two strains it was a phage/prophage transcript. This quick, easy, inexpensive, and reproducible method will facilitate the presentation of operons, transcripts, and UTR predictions alongside CDS and protein predictions in bacterial genome annotation as important resources for the research community.

Common analysis of direct RNA sequencinG CUrrently leads to misidentification of m5C at GCU motifs.

RNA modifications, such as methylation, can be detected with Oxford Nanopore Technologies direct RNA sequencing. One commonly used tool for detecting 5-methylcytosine (m5C) modifications is Tombo, which uses an "Alternative Model" to detect putative modifications from a single sample. We examined direct RNA sequencing data from diverse taxa including viruses, bacteria, fungi, and animals. The algorithm consistently identified a m5C at the central position of a GCU motif. However, it also identified a m5C in the same motif in fully unmodified in vitro transcribed RNA, suggesting that this is a frequent false prediction. In the absence of further validation, several published predictions of m5C in a GCU context should be reconsidered, including those from human coronavirus and human cerebral organoid samples.

Immunogenomic profile at baseline predicts host susceptibility to clinical malaria.

Introduction: Host gene and protein expression impact susceptibility to clinical malaria, but the balance of immune cell populations, cytokines and genes that contributes to protection, remains incompletely understood. Little is known about the determinants of host susceptibility to clinical malaria at a time when acquired immunity is developing. Methods: We analyzed peripheral blood mononuclear cells (PBMCs) collected from children who differed in susceptibility to clinical malaria, all from a small town in Mali. PBMCs were collected from children aged 4-6 years at the start, peak and end of the malaria season. We characterized the immune cell composition and cytokine secretion for a subset of 20 children per timepoint (10 children with no symptomatic malaria age-matched to 10 children with >2 symptomatic malarial illnesses), and gene expression patterns for six children (three per cohort) per timepoint. Results: We observed differences between the two groups of children in the expression of genes related to cell death and inflammation; in particular, inflammatory genes such as CXCL10 and STAT1 and apoptotic genes such as XAF1 were upregulated in susceptible children before the transmission season began. We also noted higher frequency of HLA-DR+ CD4 T cells in protected children during the peak of the malaria season and comparable levels cytokine secretion after stimulation with malaria schizonts across all three time points. Conclusion: This study highlights the importance of baseline immune signatures in determining disease outcome. Our data suggests that differences in apoptotic and inflammatory gene expression patterns can serve as predictive markers of susceptibility to clinical malaria.

Common Analysis of Direct RNA SequencinG CUrrently Leads to Misidentification of 5-Methylcytosine Modifications at GCU Motifs.

RNA modifications, such as méthylation, can be detected with Oxford Nanopore Technologies direct RNA sequencing. One commonly used tool for detecting 5-methylcytosine (m5C) modifications is Tombo, which uses an "Alternative Model" to detect putative modifications from a single sample. We examined direct RNA sequencing data from diverse taxa including virus, bacteria, fungi, and animals. The algorithm consistently identified a 5-methylcytosine at the central position of a GCU motif. However, it also identified a 5-methylcytosine in the same motif in fully unmodified in vitro transcribed RNA, suggesting that this a frequent false prediction. In the absence of further validation, several published predictions of 5-methylcytosine in human coronavirus and human cerebral organoid RNA in a GCU context should be reconsidered.

Accumulation of endosymbiont genomes in an insect autosome followed by endosymbiont replacement.

Eukaryotic genomes can acquire bacterial DNA via lateral gene transfer (LGT).1 A prominent source of LGT is Wolbachia,2 a widespread endosymbiont of arthropods and nematodes that is transmitted maternally through female germline cells.3,4 The DNA transfer from the Wolbachia endosymbiont wAna to Drosophila ananassae is extensive5-7 and has been localized to chromosome 4, contributing to chromosome expansion in this lineage.6 As has happened frequently with claims of bacteria-to-eukaryote LGT, the contribution of wAna transfers to the expanded size of D. ananassae chromosome 4 has been specifically contested8 owing to an assembly where Wolbachia sequences were classified as contaminants and removed.9 Here, long-read sequencing with DNA from a Wolbachia-cured line enabled assembly of 4.9 Mbp of nuclear Wolbachia transfers (nuwts) in D. ananassae and a 24-kbp nuclear mitochondrial transfer. The nuwts are <8,000 years old in at least two locations in chromosome 4 with at least one whole-genome integration followed by rapid extensive duplication of most of the genome with regions that have up to 10 copies. The genes in nuwts are accumulating small indels and mobile element insertions. Among the highly duplicated genes are cifA and cifB, two genes associated with Wolbachia-mediated Drosophila cytoplasmic incompatibility. The wAna strain that was the source of nuwts was subsequently replaced by a different wAna endosymbiont. Direct RNA Nanopore sequencing of Wolbachia-cured lines identified nuwt transcripts, including spliced transcripts, but functionality, if any, remains elusive.

X-treme loss of sequence diversity linked to neo-X chromosomes in filarial nematodes.

The sequence diversity of natural and laboratory populations of Brugia pahangi and Brugia malayi was assessed with Illumina resequencing followed by mapping in order to identify single nucleotide variants and insertions/deletions. In natural and laboratory Brugia populations, there is a lack of sequence diversity on chromosome X relative to the autosomes (πX/πA = 0.2), which is lower than the expected (πX/πA = 0.75). A reduction in diversity is also observed in other filarial nematodes with neo-X chromosome fusions in the genera Onchocerca and Wuchereria, but not those without neo-X chromosome fusions in the genera Loa and Dirofilaria. In the species with neo-X chromosome fusions, chromosome X is abnormally large, containing a third of the genetic material such that a sizable portion of the genome is lacking sequence diversity. Such profound differences in genetic diversity can be consequential, having been associated with drug resistance and adaptability, with the potential to affect filarial eradication.

Pim kinase inhibitor co-treatment decreases alternative non-homologous end-joining DNA repair and genomic instability induced by topoisomerase 2 inhibitors in cells with FLT3 internal tandem duplication.

Acute myeloid leukemia (AML) with fms-like tyrosine kinase 3 internal tandem duplication (FLT3-ITD) relapses with new chromosome abnormalities following chemotherapy, implicating genomic instability. Error-prone alternative non-homologous end-joining (Alt-NHEJ) DNA double-strand break (DSB) repair is upregulated in FLT3-ITD-expresssing cells, driven by c-Myc. The serine/threonine kinase Pim-1 is upregulated downstream of FLT3-ITD, and inhibiting Pim increases topoisomerase 2 (TOP2) inhibitor chemotherapy drug induction of DNA DSBs and apoptosis. We hypothesized that Pim inhibition increases DNA DSBs by downregulating Alt-NHEJ, also decreasing genomic instability. Alt-NHEJ activity, measured with a green fluorescent reporter construct, increased in FLT3-ITD-transfected Ba/F3-ITD cells treated with TOP2 inhibitors, and this increase was abrogated by Pim kinase inhibitor AZD1208 co-treatment. TOP2 inhibitor and AZD1208 co-treatment downregulated cellular and nuclear expression of c-Myc and Alt-NHEJ repair pathway proteins DNA polymerase θ, DNA ligase 3 and XRCC1 in FLT3-ITD cell lines and AML patient blasts. ALT-NHEJ protein downregulation was preceded by c-Myc downregulation, inhibited by c-Myc overexpression and induced by c-Myc knockdown or inhibition. TOP2 inhibitor treatment increased chromosome breaks in metaphase spreads in FLT3-ITD-expressing cells, and AZD1208 co-treatment abrogated these increases. Thus Pim kinase inhibitor co-treatment both enhances TOP2 inhibitor cytotoxicity and decreases TOP2 inhibitor-induced genomic instability in cells with FLT3-ITD.

Comparison of long-read sequencing technologies in interrogating bacteria and fly genomes.

The newest generation of DNA sequencing technology is highlighted by the ability to generate sequence reads hundreds of kilobases in length. Pacific Biosciences (PacBio) and Oxford Nanopore Technologies (ONT) have pioneered competitive long read platforms, with more recent work focused on improving sequencing throughput and per-base accuracy. We used whole-genome sequencing data produced by three PacBio protocols (Sequel II CLR, Sequel II HiFi, RS II) and two ONT protocols (Rapid Sequencing and Ligation Sequencing) to compare assemblies of the bacteria Escherichia coli and the fruit fly Drosophila ananassae. In both organisms tested, Sequel II assemblies had the highest consensus accuracy, even after accounting for differences in sequencing throughput. ONT and PacBio CLR had the longest reads sequenced compared to PacBio RS II and HiFi, and genome contiguity was highest when assembling these datasets. ONT Rapid Sequencing libraries had the fewest chimeric reads in addition to superior quantification of E. coli plasmids versus ligation-based libraries. The quality of assemblies can be enhanced by adopting hybrid approaches using Illumina libraries for bacterial genome assembly or polishing eukaryotic genome assemblies, and an ONT-Illumina hybrid approach would be more cost-effective for many users. Genome-wide DNA methylation could be detected using both technologies, however ONT libraries enabled the identification of a broader range of known E. coli methyltransferase recognition motifs in addition to undocumented D. ananassae motifs. The ideal choice of long read technology may depend on several factors including the question or hypothesis under examination. No single technology outperformed others in all metrics examined.

Complete Mitochondrial Genome Sequence of Mansonella perstans.

The 13,647-bp complete mitochondrial genome of Mansonella perstans was sequenced and is syntenic to the mitochondrial genome of Mansonella ozzardi Phylogenetic analysis of the mitochondrial genome is consistent with the known phylogeny of ONC5 group filarial nematodes.

Nearly Complete Genome Sequence of Brugia pahangi FR3.

Brugia pahangi is a zoonotic parasite that is closely related to human-infecting filarial nematodes. Here, we report the nearly complete genome of Brugia pahangi, including assemblies of four autosomes and an X chromosome, with only seven gaps. The Y chromosome is still not completely assembled.

Complete Genome Sequence of wBp, the Wolbachia Endosymbiont of Brugia pahangi FR3.

Lymphatic filariasis is a devastating disease caused by filarial nematode roundworms, which contain obligate Wolbachia endosymbionts. Here, we assembled the genome of wBp, the Wolbachia endosymbiont of the filarial nematode Brugia pahangi, from Illumina, Pacific Biosciences, and Oxford Nanopore data. The complete, circular genome is 1,072,967 bp.

Drug Repurposing of Bromodomain Inhibitors as Potential Novel Therapeutic Leads for Lymphatic Filariasis Guided by Multispecies Transcriptomics.

To better understand the transcriptomic interplay of organisms associated with lymphatic filariasis, we conducted multispecies transcriptome sequencing (RNA-Seq) on the filarial nematode Brugia malayi, its Wolbachia endosymbiont wBm, and its laboratory vector Aedes aegypti across the entire B. malayi life cycle. In wBm, transcription of the noncoding 6S RNA suggests that it may be a regulator of bacterial cell growth, as its transcript levels correlate with bacterial replication rates. For A. aegypti, the transcriptional response reflects the stress that B. malayi infection exerts on the mosquito with indicators of increased energy demand. In B. malayi, expression modules associated with adult female samples consistently contained an overrepresentation of genes involved in chromatin remodeling, such as the bromodomain-containing proteins. All bromodomain-containing proteins encoded by B. malayi were observed to be upregulated in the adult female, embryo, and microfilaria life stages, including 2 members of the bromodomain and extraterminal (BET) protein family. The BET inhibitor JQ1(+), originally developed as a cancer therapeutic, caused lethality of adult worms in vitro, suggesting that it may be a potential therapeutic that can be repurposed for treating lymphatic filariasis.IMPORTANCE The current treatment regimen for lymphatic filariasis is mostly microfilaricidal. In an effort to identify new drug candidates for lymphatic filariasis, we conducted a three-way transcriptomics/systems biology study of one of the causative agents of lymphatic filariasis, Brugia malayi, its Wolbachia endosymbiont wBm, and its vector host Aedes aegypti at 16 distinct B. malayi life stages. B. malayi upregulates the expression of bromodomain-containing proteins in the adult female, embryo, and microfilaria stages. In vitro, we find that the existing cancer therapeutic JQ1(+), which is a bromodomain and extraterminal protein inhibitor, has adulticidal activity in B. malayi.

Complete Genome Sequence of wAna, the Wolbachia Endosymbiont of Drosophila ananassae.

Here, we present the complete genome sequence of the Wolbachia endosymbiont wAna, isolated from Drosophila ananassae and derived from Oxford Nanopore and Illumina sequencing. We anticipate that this will aid in Wolbachia comparative genomics and the assembly of D. ananassae specifically in regions containing extensive lateral gene transfer events.

TwinBLAST: When Two Is Better than One.

Analysis of sequence read pairs can be essential for characterizing structural variation, including junction-spanning pairs of reads (JSPRs) suggesting recent lateral/horizontal gene transfer. TwinBLAST can be used to facilitate this analysis of JSPRs by enabling the visualization and curation of two BLAST reports side by side in a single interface.

Multispecies Transcriptomics Data Set of Brugia malayi, Its Wolbachia Endosymbiont wBm, and Aedes aegypti across the B. malayi Life Cycle.

Here, we present a comprehensive transcriptomics data set of Brugia malayi, its Wolbachia endosymbiont wBm, and its vector host. This study samples from 16 stages across the entire B. malayi life cycle, including stage 1 through 4 larvae, adult males and females, embryos, immature microfilariae, and mature microfilariae.

Targeted enrichment outperforms other enrichment techniques and enables more multi-species RNA-Seq analyses.

Enrichment methodologies enable the analysis of minor members in multi-species transcriptomic data. We compared the standard enrichment of bacterial and eukaryotic mRNA to a targeted enrichment using an Agilent SureSelect (AgSS) capture for Brugia malayi, Aspergillus fumigatus, and the Wolbachia endosymbiont of B. malayi (wBm). Without introducing significant systematic bias, the AgSS quantitatively enriched samples, resulting in more reads mapping to the target organism. The AgSS-enriched libraries consistently had a positive linear correlation with their unenriched counterparts (r2 = 0.559-0.867). Up to a 2,242-fold enrichment of RNA from the target organism was obtained following a power law (r2 = 0.90), with the greatest fold enrichment achieved in samples with the largest ratio difference between the major and minor members. While using a single total library for prokaryote and eukaryote enrichment from a single RNA sample could be beneficial for samples where RNA is limiting, we observed a decrease in reads mapping to protein coding genes and an increase in multi-mapping reads to rRNAs in AgSS enrichments from eukaryotic total RNA libraries compared to eukaryotic poly(A)-enriched libraries. Our results support a recommendation of using AgSS targeted enrichment on poly(A)-enriched libraries for eukaryotic captures, and total RNA libraries for prokaryotic captures, to increase the robustness of multi-species transcriptomic studies.

Lateral gene transfer between prokaryotes and eukaryotes.

Lateral gene transfer (LGT) is an all-encompassing term for the movement of DNA between diverse organisms. LGT is synonymous with horizontal gene transfer, and the terms are used interchangeably throughout the scientific literature. While LGT has been recognized within the bacteria domain of life for decades, inter-domain LGTs are being increasingly described. LGTs between bacteria and complex multicellular organisms are of interest because they challenge the long-held dogma that such transfers could only occur in closely-related, single-celled organisms. Scientists will continue to challenge our understanding of LGT as we sequence more, diverse organisms, as we sequence more endosymbiont-colonized arthropods, and as we continue to appreciate LGT events, both young and old.

Extensive duplication of the Wolbachia DNA in chromosome four of Drosophila ananassae.

BACKGROUND: Lateral gene transfer (LGT) from bacterial Wolbachia endosymbionts has been detected in ~20% of arthropod and nematode genome sequencing projects. Many of these transfers are large and contain a substantial part of the Wolbachia genome. RESULTS: Here, we re-sequenced three D. ananassae genomes from Asia and the Pacific that contain large LGTs from Wolbachia. We find that multiple copies of the Wolbachia genome are transferred to the Drosophila nuclear genome in all three lines. In the D. ananassae line from Indonesia, the copies of Wolbachia DNA in the nuclear genome are nearly identical in size and sequence yielding an even coverage of mapped reads over the Wolbachia genome. In contrast, the D. ananassae lines from Hawaii and India show an uneven coverage of mapped reads over the Wolbachia genome suggesting that different parts of these LGTs are present in different copy numbers. In the Hawaii line, we find that this LGT is underrepresented in third instar larvae indicative of being heterochromatic. Fluorescence in situ hybridization of mitotic chromosomes confirms that the LGT in the Hawaii line is heterochromatic and represents ~20% of the sequence on chromosome 4 (dot chromosome, Muller element F). CONCLUSIONS: This collection of related lines contain large lateral gene transfers composed of multiple Wolbachia genomes that constitute >2% of the D. ananassae genome (~5 Mbp) and partially explain the abnormally large size of chromosome 4 in D. ananassae.

Large-scale intersubspecific recombination in the plant-pathogenic bacterium Xylella fastidiosa is associated with the host shift to mulberry.

Homologous recombination plays an important role in the structuring of genetic variation of many bacteria; however, its importance in adaptive evolution is not well established. We investigated the association of intersubspecific homologous recombination (IHR) with the shift to a novel host (mulberry) by the plant-pathogenic bacterium Xylella fastidiosa. Mulberry leaf scorch was identified about 25 years ago in native red mulberry in the eastern United States and has spread to introduced white mulberry in California. Comparing a sequence of 8 genes (4,706 bp) from 21 mulberry-type isolates to published data (352 isolates representing all subspecies), we confirmed previous indications that the mulberry isolates define a group distinct from the 4 subspecies, and we propose naming the taxon X. fastidiosa subsp. morus. The ancestry of its gene sequences was mixed, with 4 derived from X. fastidiosa subsp. fastidiosa (introduced from Central America), 3 from X. fastidiosa subsp. multiplex (considered native to the United States), and 1 chimeric, demonstrating that this group originated by large-scale IHR. The very low within-type genetic variation (0.08% site polymorphism), plus the apparent inability of native X. fastidiosa subsp. multiplex to infect mulberry, suggests that this host shift was achieved after strong selection acted on genetic variants created by IHR. Sequence data indicate that a single ancestral IHR event gave rise not only to X. fastidiosa subsp. morus but also to the X. fastidiosa subsp. multiplex recombinant group which infects several hosts but is the only type naturally infecting blueberry, thus implicating this IHR in the invasion of at least two novel native hosts, mulberry and blueberry.

Recent evolutionary radiation and host plant specialization in the Xylella fastidiosa subspecies native to the United States.

The bacterial pathogen, Xylella fastidiosa, infects many plant species in the Americas, making it a good model for investigating the genetics of host adaptation. We used multilocus sequence typing (MLST) to identify isolates of the native U.S. subsp. multiplex that were largely unaffected by intersubspecific homologous recombination (IHR) and to investigate how their evolutionary history influences plant host specialization. We identified 110 "non-IHR" isolates, 2 minimally recombinant "intermediate" ones (including the subspecific type), and 31 with extensive IHR. The non-IHR and intermediate isolates defined 23 sequence types (STs) which we used to identify 22 plant hosts (73% trees) characteristic of the subspecies. Except for almond, subsp. multiplex showed no host overlap with the introduced subspecies (subspecies fastidiosa and sandyi). MLST sequences revealed that subsp. multiplex underwent recent radiation (<25% of subspecies age) which included only limited intrasubspecific recombination (ρ/θ = 0.02); only one isolated lineage (ST50 from ash) was older. A total of 20 of the STs grouped into three loose phylogenetic clusters distinguished by nonoverlapping hosts (excepting purple leaf plum): "almond," "peach," and "oak" types. These host differences were not geographical, since all three types also occurred in California. ST designation was a good indicator of host specialization. ST09, widespread in the southeastern United States, only infected oak species, and all peach isolates were ST10 (from California, Florida, and Georgia). Only ST23 had a broad host range. Hosts of related genotypes were sometimes related, but often host groupings crossed plant family or even order, suggesting that phylogenetically plastic features of hosts affect bacterial pathogenicity.